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\  ABSTRACT 

The  instability  of  hypersonic  boundary-layer  flows  over  flat  plates  is  considered.  The 
viscosity  of  the  fluid  is  taken  to  be  governed  by  Sutherland’s  law,  which  gives  a  much  more 
accurate  representation  of  the  temperature  dependence  of  fluid  viscosity  at  hypersonic  speeds 
than  Chapman’s  approximate  linear  law;  although  at  lower  speeds  the  temperature  variation 
of  the  mean  state  is  less  pronounced  so  that  the  Chapman  law  can  be  used  with  some  con¬ 
fidence.  Attention  is  focussed  on  the  so-called  "^orticity^^ode  of  instability  of  the  viscous 
hypersonic  boundary  layer.-.  This  is  thought  to  be  the  fastest  growing  inviscid  disturbance 
at  hypersonic  speeds;  it  is  tfiso  believed  to  have  an  asymptotically  larger  growth  rate  than 
any  viscous  or  centrifugal  insfaHIG^.  As  a  starting  point  we  investigat^he  instability  of 
the  hypersonic  boundary  layer  exists  far  downstream  from  the  leading  edge  of  the 

plate.  In  this  regime  the  shock  that  is  attached  to  the  leading  edge  of  the  plate  plays  no 
role,  so  that  the  basic  boundary  layer  is  non-interactive,  .ft*  is  Huuwu  lln^  the  vortidty  mode 
of  instability  of  this  flow  operates  on  a  significantly  different  lengthscale  than  that  obtained 
if  a  Chapman  viscosity  law  is  assume^-see-Snut^h  and  Brown,  1989)*.  In  particular,  we  find 
that^he  growth  rate  predicted  by  a  linear  viscosity  law  overestimates  the  size  of  the  growth 
rate  by  0(^*1.  Next,  the  development  of  the  vortidty  mode  as  the  wavenumber  decreases  is 
described,  ana  it  is  shown  that  acoustic  modes  emerge  when  the  wavenumber  has  decreased 

from  it’s  Ofl)  initial  value  to  Finally,  the  inviscid  instability  of  the  boundary 

layer  near  tne  leading  edge  in  the  intera^ion  zone  is  discussed  and  particular  attention  is 
focussed  on  the  strong  interaction  region  v^hich  occurs  sufficiently  close  to  the  leading  edge. 
We  find  that  the  vorticity  mode  in  this  regjime  is  again  unstable,  and  that  it  is  concentrated 
in  the  transition  layer  at  the  edge  of  the  boundary  layer  where  the  temperature  adjusts 
from  its  large,  O(M^),  value  in  the  viscous  boundary  layer,  to  its  0(1)  free  stream  value. 
The  existence  of  the  shock  indirectly,  but  significantly,  influences  the  instability  problem  by 
modifying  the  basic  flow  structure  in  this  layer. 

^  r  I 

2''  ^  ■  I  ..t  -- 
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1.  INTRODUCTION 

Our  concern  is  with  the  inviscid  mode  of  instability  of  hypersonic  boundary-layer  flows 
over  flat  plates.  In  the  first  instance  we  will  consider  a  regime  where  there  are  no  shocks 
present,  and  then  we  will  show  how  the  instability  problem  is  significantly  modified  by 
the  presence  of  a  shock  in  the  flow  field.  The  motivation  for  this  and  related  work  on 
hypersonic  boundary-layer  instability  theory  is  the  renewed  interest  in  hypersonic  flight 
which  has  been  stimulated  by  plans  to  build  a  successor  to  the  Space  Shuttle.  A  primary 
concern  with  such  a  vehicle  is  the  question  of  where  transition  will  occur  over  a  wide 
range  of  Mach  numbers  and  whether  it  can  be  controlled.  At  the  largest  relevant  Mach 
numbers,  say  Mach  20-25,  the  extremely  high  temperatures  associated  with  the  flow  would 
destroy  the  vehicle  unless  it  were  cooled,  so  that  it  is  of  interest  to  know  the  effect  of 
the  wall  temperature  on  the  instability  properties  of  the  flow.  The  purpose  of  this  paper 
is  to  determine  the  inviscid  instability  characteristics  of  physically  realistic  hypersonic 
boundary-layer  flows.  We  note  here  in  passing  that  there  is  a  simple  generalization  of 
Rayleigh’s  (incompressible)  inflection  point  theorem  to  compressible  flows  (Lees  &  Lin, 
1946),  and  that  many  compressible  boundary  layers  turn  out  to  be  inviscidly  unstable  even 
though  their  incompressible  counterparts  are  stable.  This  is  a  significant  result  because 
the  growth  rates  of  inviscid  disturbances  are  generally  much  larger  than  those  of  viscous 
or  centrifugal  instabilities  so  that  they  are  the  likely  cause  of  transition  to  turbulence 
in  most  situations.  The  modes  which  we  discuss  in  this  paper  can  be  referred  to  as 
generalized  inflection  point  modes  because,  when  neutral,  their  phase  speed  is  equal  to  the 
fluid  velocity  at  the  generalized  inflection  point.  Furthermore  the  eigenfunctions  of  the 
modes  are  localized  around  that  point. 

For  convenience  we  will  concentrate  on  high-Reynolds-number  flow  past  a  flat  plate, 
although  many  aspects  of  our  analysis  are  applicable  to  other  boundary-layer  flows  (e.g. 
flow  past  a  wedge).  The  underlying  steady  flows  that  we  study  depend  upon  the  Mach 
number,  the  Prandtl  number  and  the  choice  of  the  viscosity  law;  the  complications  arising 
from  real  gas  effects  are  not  investigated.  Throughout  we  assume  that  the  fluid  viscosity 
is  Newtonian  and  is  adequately  described  by  Sutherland’s  formula.  We  will  also  take 
the  Prandtl  number  to  be  one,  noting  that  it  is  relatively  straightforward  to  relax  this 
restriction.  Moreover  the  relaxation  of  the  latter  restriction  does  not  significantly  alter  the 
qusditative  features  of  the  results  we  present  here. 
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Reshotko  (1976)  and  Mack  (1987)  have  reviewed  earlier  work  on  the  linear  instability 
of  high-Reynolds-number  compressible  flows.  Many  of  these  studies  are  based  on  the  Orr- 
Sommerfeld  equation;  for  a  critique  of  the  mathematical  rigor  of  this  approach  see  Smith 
(1979, 1989).  Here  we  examine  the  linear  stability  of  high-Reynolds-number  flows  by  means 
of  formal  asymptotic  expansions;  for  example  Smith  (1989)  has  applied  triple-deck  theory 
to  the  lower-branch  viscous  Tollmien-Schlichting  modes  of  compressible  boundary  layers. 
Seddougui,  Bowles  &  Smith  (1989)  have  extended  this  theory  to  include  the  effects  of 
severe  wall  cooling,  while  Cowley  &  Hall  (1988)  have  shown  how  such  modes  can  interact 
with  a  shock  at  large  Mach  number.  However,  the  viscous  modes  have  relatively  small 
growth  rates,  and  our  main  concern  here  will  be  with  inviscid  modes.  The  nature  of  the 
asymptotic  expansion  procedure  in  these  investigations  clearly  depends  on  the  nature  of 
the  mode  of  instability.  In  fact  the  third  type  of  instability  responsible  for  boundary- 
layer  transition,  the  Gortler  vortex  mode,  develops  an  asymptotic  structure  at  high  Mach 
numbers  closely  related  to  that  found  here;  see  Hall  and  Fu  (1989). 

When  a  quasi-parallel  approximation  is  formally  justifiable  because  the  Reynolds  num¬ 
ber  is  large,  inviscid  modes  satisfy  the  compressible  Rayleigh  equation.  Numerical  solutions 
to  this  equation  have  been  reported  by,  inter  alia,  Mack  (1984,  1987)  for  boundary-layer 
flows,  Grosch  &  Jackson  (1989)  for  shear  flows,  and  Papageorgiou  (1989)  for  wake  flows. 
For  fluids  satisfying  a  Chapman  viscosity  law,  high-Mach-number  asymptotic  solutions  to 
this  equation  for  the  so-called  “acoustic”  boundary-layer  modes  have  been  obtained  by 
Cowley  Sc  Hall  (1989),  while  Smith  Sc  Brown  (1989)  have  identified  the  asymptotic  form 
of  the  “vorticity”  mode.  Goldstein  Sc  Balsa  (1989)  have  given  an  asymptotic  solution  for 
high-Mach-number,  shear-layer  modes  of  instability  for  a  Chapman  fluid.  Though  the  basic 
states  investigated  by  the  two  latter  pairs  of  authors  are  different,  they  found  essentailly 
the  same  most  unstable  eigenvalue  because  it  corresponds  to  a  disturbance  trapped  in  a 
thin  layer  where  the  overall  features  of  the  basic  state  are  unimportant.  However  Goldstein 
and  Balsa  did  not  spot  the  exact  solution  of  the  vorticity  mode  equation  found  by  Smith 
and  Brown. 

Both  in  the  above  mentioned  boundary-layer  analyses,  and  the  hypersonic  Gortler 
vortex  instability  analysis  of  Hall  Sc  Fu  (1989),  one  of  the  key  asymptotic  regions  for  the 
case  of  a  Chapman  viscosity  law  is  a  logarithmically  thin  layer  which  develops  due  to  the 
exponential  decay  of  the  underlying  steady  temperature  field  away  from  the  wall.  However, 
Chapman’s  viscosity  law  is  not  exact,  and  was  introduced  as  a  useful  interpolation  law 
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which  greatly  simplified  steady  boundary-layer  calculations  (e.g.  see  Stewartson  1964).  At 
the  large  temperatures  typical  in  hypersonic  flows,  it  differs  significantly  from  the  more 
precise  Sutherland’s  formula.  In  fact  Chapman’s  law  is  simply  a  linear  approximation  to 
the  viscosity-temperature  dependence  of  the  fluid;  it  is  therefore  of  questionable  validity 
in  the  hypersonic  limit.  At  high  Mach  numbers  the  steady  temperature  field  in  a  fluid 
satisfying  Sutherland’s  formula  initially  decays  algebraically  away  from  the  wall,  before 
reverting  to  exponential  decay  in  an  asymptotic  region  “far”  from  the  wall  (e.g.  Freeman 
&  Lam  1959).  The  effect  of  this  algebraic  decay  is  to  change  significantly  the  scalings  in  the 
transition  region;  in  particular  the  asymptotic  expansions  proceed  in  inverse  powers  of  M 
rather  than  •\/log(M).  Moreover  the  wavelength  of  the  most  unstable  disturbance  varies 
by  a  factor  of  •v/log(M)  in  the  two  cases.  We  note  that  a  similar  difference  in  scalings  is 
evident  in  the  interaction  region  of  steady  hypersonic  flow  past  a  flat  plate.  In  that  case 
Lee  &  Cheng  (1969)  have  shown  that  the  shock-heating  transition  layer  is  logarithmically 
thin  for  Chapman’s  viscosity  law,  whereas  for  a  power-law  viscosity  formula,  and  hence 
for  Sutherland’s  formula,  the  scaling  for  the  transition  layer  is  algebraic  (Bush  1966). 

The  flows  which  we  consider  here  are  appropriate  to  different  stages  of  hypersonic  flow 
past  a  semi-infinite  flat  plate.  In  the  first  instance  we  shall  consider  the  instability  of  a 
non-interactive  flow  which  is  appropriate  to  large  distances  downstream  of  the  leading  edge 
of  the  plate.  Here  the  attached  shock  at  the  leading  edge  has  no  effect  on  the  flow  field  and 
the  basic  state  is  the  Sutherland  law  counterpart  of  that  discussed  by  Smith  and  Brown 
(1989).  This  basic  state,  and  the  Rayleigh  equation  which  governs  the  inviscid  instability, 
are  discussed  in  Section  2.  The  dispersion  relationship  associated  with  this  equation  is 
then  solved  in  Section  3.  We  shall  discuss  the  growth  rate  of  the  mode  over  the  whole 
range  of  unstable  wavenumbers  and  discuss  how  the  disturbance  is  related  to  the  acoustic 
mode  at  small  wavenumbers. 

Then  in  Section  4  we  go  on  to  discuss  the  basic  state  in  the  so-called  “strong  interaction 
region”  further  upstream.  The  description  of  the  basic  flow  in  this  regime  for  a  power  law 
fluid  is  due  to  Luniev  (1959)  and  Bush  (1966).  We  shall  discuss  how  the  interactive  system 
for  the  flow  can  be  formulated  when  the  Mach  number  is  of  order  iZK  However  this 
system  can  only  be  solved  numerically  and  we  are  not  aware  of  any  published  results  on 
that  problem.  Nevertheless  we  can  of  course  still  consider  the  instability  of  that  state, 
and  we  derive  the  appropriate  (quasi-parallel)  stability  equations  in  this  regime.  The 
strong  hypersonic  interaction  limit  then  corresponds  to  taking  the  further  limit  of  the 
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streamwise  variable  tending  to  zero.  In  that  case  a  similarity  solution  for  the  basic  flow 
can  be  found,  Bush  (1966),  and  a  re-scaled  Rayleigh  equation  for  the  disturbance  is  found. 
The  solution  of  that  equation  is  discussed  in'  Section  5.  We  could  have  instead  considered 
the  weak  hypersonic  limit  further  downstream  where  Bush  and  Cross  (1967)  have  given  an 
appropriate  asymptotic  description.  We  choose  to  concentrate  on  the  strong  interaction 
regime  because  it  is,  to  a  certain  extent,  simpler.  Further,  if  the  flow  is  unstable  in  this 
regime  it  is  likely  that  growing  disturbances  will  originate  here.  Finally  in  Section  6  we 
shall  draw  some  conclusions. 
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2.  NON-INTERACTIVE  STEADY  FLOWS 

We  begin  by  considering  the  stability  of  steady  hypersonic  flow  far  downstream  from 
any  leading-edge  interaction  region;  in  particular,  if  L  is  the  distance  from  the  leading 
edge,  and  U^o  >  doo »  Poo  and  ,  are  the  velocity,  sound  speed,  density  and  shear  viscosity 
of  the  free  stream  flow,  then  we  assume  that  the  Reynolds  number. 


Po.U^L 


(2.1a) 


is  larger  than  whatever  power  of  the  Mach  number. 


M  = 


(2.16) 


is  necessary  for  interactive  and/or  non-parallel  effects  to  be  negligible  (see  below  for  a  more 
precise  restriction).  We  adopt  a  non-dimensionalisation  based  on  coordinates  Lx  (where  x 
is  in  the  direction  of  flow  and  y  is  normal  to  plate),  velocities  UooM.,  time  Lt/Uoo  »  pressure 
PoaU^p,  density  PooP>  temperature  TooT,  and  shear  and  bulk  viscosities  PooP  and  PooP^ 
respectively,  where  the  subscript  oo  denotes  the  value  of  the  quantity  in  the  free-stream. 
On  the  assumption  that  the  fluid  is  a  perfect  gas  with  a  constant  ratio  of  specific  heats  7, 
the  governing  equations  of  the  flow  are 


+  V  ■  (pu)  =  0, 

p^  =  -Sp+  |[2V  •  (/<£)  +  -  5m)V  ■  Jt)|, 


DT  .  \%r2^P  f 


oDp  1  (7  —  1)M'' 

■  (.pST)  +  „ 


where 


7M>  =  pT, 


2  ^  dxj  dxi '  ’ 


$  =  2p| :  I  -I-  (p'  -  2P)(  V  ■  u)^. 


(2.2a) 


(2.26) 


(2.2c) 

(2.2<f) 


(2.3a) 


(2.36) 


and  Pr  is  the  constant  Prandtl  number.  For  a  shear  viscosity  obeying  Sutherland’s  law 


=  fi±£)  Ti 

Vr  +  c/  ’ 


{2,3c) 
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where  C 


110.4 


for  air  temperatures  measured  in  degrees  Kelvin.  In  the  numerical  cal¬ 


culations  discussed  in  the  next  section  we  took  too  =  216.9.  The  boundary-layer  equation 
can  be  recovered  by  first  substituting 


R-^V, 


(2.4a,  6) 


where  the  Dorodnitsyn-Howarth  variable,  C,  is  introduced  for  convenience,  and  then  taking 
the  limit  R  —*  oo. 

For  steady  two-dimensional  flow  over  a  fiat  plate,  a  similarity  solution  to  these  equa¬ 
tions  exists.  With 


^  =  ~r7T==^'  u  =  V’(,  = -(V’.  + 

V  (1  +  C)x 

V)  =  \/(l  +  C)x  fir)),  T  =  Tit)),  p  =  pii)), 
the  governing  similarity  equations  are  found  to  be 


(2.5a) 

(2.5fe) 


pr=  1, 


subject  to  the  boundary  conditions 


/(O)  =  no)  =  0,  7'(oo)  =  r(oo)  =  1, 


(2.6a) 

(2.66) 

(2.6c) 


(2.6d) 


and  r(0)  =  r„  (fixed  wall  —  temperature),  or  r'(0)  =  0  (insulated  wall).  (2.6e) 

For  simplicity  we  will  focus  attention  on  Pr  =  1,  and  denote  by  T,  the  wall  temperature 
when  the  boundary  is  insulated.  Then,  as  is  well  known  (e.g.  Stewartson  1964),  the  energy 
equation  can  be  integrated  to  yield 

r  =  1  +  ({T,  - 1)  + 1(7  -  +  /'))(!  -  /'),  (2.7) 

where  r„,  =  TtT,  and  T,  =  1  -f-  ^(7  —  1)M^. 
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The  solution  to  (2.6)  in  the  limit  of  large  Mach  number  has  been  examined  by  Freeman 
&  Lam  (1959).  They  showed  that  two  asymptotic  regions  develop,  defined  as  where  the 
coordinates  rj  and  (  =  M^rj,  respectively,  are  order  one. 


e  =  Q(i) 

In  this  region  we  write  /  =  M“>/o(0  +  •  •  •?  then  using  (2.7)  it  follows  from  (2.6b) 


+  (:7^)  !_/,)*)  =  “■ 


(2.8a) 


/o(0)  =  /6(0)  =  0, 


(2.86) 


24  723 

/o  ~  f  +  ^  1)(7  -  l)e  ~  (T.  +  1)(7  -  Ije-l  +  ■  ■  •  “  f 

We  note  that  as  ^  »  cx),  then  fo  —*  1  algebraically.  This  is  of  course  different  than 

the  corresponding  result  for  a  Chapman  fluid,  in  that  case  the  correction  term  to  the 
free-stream  speed  is  exponentially  small.  This  difierence  is  significant  because  it  leads  to 
inviscid  instabilities  which  are  unstable  at  very  different  wavelengths. 


V  =  Q(i) 

Here,  we  write 

then  gi  satisfies  the  equation 


/  =  ^  +  +  •  •  •  , 


„,J  n/i - 1(7 - I)(r>  +  ' 

\  1  +  C  -  ^{-r  -  i)(T.  +  l)9i 


(2.10a) 


subject  to 


(rt  +  l)(7  -  1)773  91-^0  as  77->oo. 


(2.106) 


In  Figure  (2.1a)  we  show  the  function  /q  for  the  adiabatic  case  Tj  =  1  with  7  =  1.4,  we 
note  the  algebraic  approach  of  the  function  to  the  free-stream  speed.  In  Figure  (2.1b) 
we  show  the  transition  layer  function  gi  corresponding  to  the  same  case.  We  note  the 
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exponential  decay  of  this  function  for  large  rj.  It  is  also  worth  pointing  out  that  for  a 
Chapman  fluid  the  transition-layer  equation  corresponding  to  (2,10a)  is  linear,  and  it’s 
solution  can  be  expressed  in  terms  of  the  exponential  function.  It  was  this  simpliflcation 
that  enabled  Smith  and  Brown  (1989)  to  spot  the  exact  solution  of  the  neutral  vorticity 
mode  in  their  study  of  the  instability  problem  in  this  layer.  We  shall  see  below  that  for 
Sutherland’s  law  we  are  not  able  to  And  a  similar  exact  solution  of  the  stability  equation. 

Sufficiently  far  downstream  the  quasi-parallel  assumption  is  valid  for  inviscid  instabil¬ 
ity  modes.  It  is  then  appropriate  to  seek  perturbations  of  the  form 

u  =  fin)  +  ■■■  +  ^uiT})exp(^i\J +  Pz  -  ct)j  +  ...,  (2.11a) 


P  = 


7M2 


-I- ...  -I-  Ap(T?)exp  I  i\ 


R 


{ax  +  Pz  -  ct)j  +  . . . , 


(2.116) 


il  +  C) 

with  similar  expressions  for  the  other  flow  quantites.  If  the  disturbance  amplitude,  A,  is 
sufficiently  small  the  pressure  perturbation  p  satisfles  the  compressible  Rayleigh  equation 


d^p 

drf 


V"  dp  „ 

-{a  +P  )T{T - — )P  =  0, 


f  -cdri 


(a2  -t-  /?2) 


(2.12a) 


where  for  convenience  we  have  rescaled  the  quantities  a,P,  c  by  dividing  them  by  x».  The 
conditions  that  there  is  no  normal  velocity  at  the  wall,  and  that  the  disturbance  is  confined 
to  the  boundary  layer,  can  be  expressed  as 


p'  =  0  on  7/  =  0,  p  0  as  7y  — ♦  oo. 


(2.126) 


Equation  (2.12a)  and  boundary  conditions  (2.126)  specify  a  temporal-stability  eigenrelation 
c  =  c{a,P)\  alternatively  from  a  spatial  stability  standpoint,  the  eigenrelation  can  be 
regarded  as  a  =  a{c,P). 
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3.  THE  FAR  DOWNSTREAM  BEHAVIOUR  OF  THE  INVISCID  MODES 

In  this  section  we  discuss  the  asymptotic  form  of  unstable  solutions  to  (2.12a)  for  the 
region  far  downstream  of  the  leading  edge  of  the  plate.  In  a  previous  investigation  Cowley 
and  Hall  (1988),  hereafter  referred  to  as  CH,  studied  the  so-called  acoustic  modes  of  (2.12a) 
in  this  region  on  the  assumption  that  the  viscosity  satisfies  Chapman’s  law.  Simultaneously 
Smith  and  Brown  (1989),  hereafter  referred  to  as  SB,  investigated  the  vorticity  mode  for 
Chapman’s  law.  The  main  difference  between  these  two  types  of  modes  is  that  the  acoustic 
modes  have  a  ~  M“^,  whilst  the  vorticity  mode,  at  least  while  it  is  close  to  neutral,  has 
a  ~  ^/ 2  log  M^.  Moreover  the  vorticity  mode  is  centred  at  the  adjustment  layer  at  the 
edge  of  the  boundary  layer  whilst  the  acoustic  one  is  concentrated  in  the  main  part  of  the 
boundary  layer. 

However,  as  indicated  above,  at  high  Mach  numbers  the  temperature  variations  in  the 
boundary  layer  are  large,  so  that  a  linear  temperature-viscosity  law  is  a  bad  approximation; 
Sutherland’s  law  should  be  used  to  give  a  better  representation  of  the  viscosity.  It  is  then 
important  to  see  how  the  asymptotic  structures  developed  by  CH  and  SB  change.  We  shall 
see  that  there  are  significant  differences. 

In  the  first  instance  we  derive  the  asymptotic  structure  of  the  solution  for  a  two- 
dimensional  vorticity  mode  of  (2.12a).  We  determine  the  neutral  values  of  a  and  c  for 
this  mode,  and  find  the  limiting  form  of  the  mode  when  the  further  limit  a  — »  0  is  taken. 
This  limiting  solution  points  to  a  sequence  of  distinguished  asymptotic  limits.  Within  the 
sequence  of  asymptotic  limits  mentioned  above  the  scaling,  a  ~  appropriate  to 

an  acoustic  mode  emerges;  we  therefore  discuss  the  latter  mode  as  a  limiting  case  of  the 
vorticity  mode. 


Q  =  0(1) 

Consider  then  the  solution  of  (2.12a)  which  has  the  eigenfunction  trapped  in  the 
temperature  adjustment  layer  at  the  edge  of  the  boundary  layer.  We  seek  a  solution  which 
has  p  =  0  and  a  =  0(1),  so  that  the  wavelength  of  the  vorticity  mode  is  comparable  with 
the  depth  of  the  adjustment  layer  in  terms  of  the  Dorodnitsyn-Howarth  variable.  From 
(2.7)  and  (2.9),  the  velocity  field,  tl,  and  temperature  field,  T,  of  the  underlying  steady 
flow  expand  as 

2G(f7) 

T=1-G  +  ..., 
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where  rj  =  and  the  function  G  =  +  1)(7  —  l)gf^  satisfies  a  rescaled  v  ersion  of 

(2.10a),  and  has  the  asymptotic  behaviour 


/j4  +  ^3-^/7  + 

(3.2a) 

G  0  exponentially  as  ^  »  00. 

(3.25) 

Here  B  is  a  constant  to  be  calculated  numerically. 
Next  we  expand  a,  c  and  p  in  the  forms 


a=^a+..., 

e=i  + _ ^ 

^  (7  -  1)(T,  +  1)M2 

p=  p+ 


C  + 


(3.3a, 6, c) 


where  we  have  assumed  that  the  disturbance  moves  downstream  with  the  fluid  speed  in 
the  adjustment  layer.  On  substituting  for  u  and  T  from  (3.1),  and  using  (3.3),  we  find 
that  the  zeroth  order  approximation  to  (2.12a)  in  the  adjustment  layer  is, 


d^p 

df\^ 


2G'  dp 
G  —  cdf) 


-  d2(l  -  Gfp  =  0. 


(3.4) 


We  note  that  three-dimensional  disturbances  satisfy  (3.4)  but  with  d  replaced  by 
d^  +  ,3^  so  that  in  the  present  regime  it  is  suflScient  for  us  to  consider  only  two-dimensional 
modes.  Equation  (3.4)  is  to  be  solved  subject  to  p  vanishing  in  the  limits  r?  — »  0  and 
Tj  oo,  so  that  the  disturbance  is  confined  to  the  adjustment  layer.  For  T)  >>  1  it  follows 
from  (3.4)  that  p  decays  like  exp(— d^),  whilst  for  7  <<  1  a  WKB  solution  of  (3.4)  can  be 
expressed  in  the  form 

p^exp[-J  0(r?)]d»7,  (3.5) 

where 

^  9d  2  Hd  , 

® - ^  +  jj3_v/7  +2Ca-|-  ...  as  fj ->  0.  (3.6) 

First,  we  restrict  our  attention  to  the  neutral  case,  c  is  then  real  and  can  be  evaluated 
by  finding  the  fluid  speed  correct  to  order  at  the  generalized  inflection  point  which 
is  located  where 

^_H5L-n 

Q,  T 
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i.e.  where 


G"  2G' 

l-G  ~ 

A  numerical  solution  to  (2.10a)  using  a  Runge-Kutta  method  shows  that  this  occurs  when 
Tj  ~  1.971510,  in  which  case 

c  =  -.993934  (3.7) 

The  corresponding  real  value  of  a  is  obtained  by  integrating  (3.4)  from  /7  =  0  to  ?7  =  oo  with 
an  appropriate  treatment  at  the  generalised  inflection  point.  Such  a  calculation  predicts 
that  the  neutral  value  of  a  is 

a  0.645065.  (3.8) 

However,  of  greater  significance  are  the  unstable  eigenrnodes.  Figure  (3.1)  illustrates 
the  dependence  of  the  growth  rate,  aCi,  on  the  real  wavenumber  a  .  We  see  that  the  growth 
rate  attains  it’s  maximum  value  of  ~  0.256853  at  d  ~  0.143619;  further  it  turns  out  that 
the  acoustic  modes  have  smaller  growth  rates  (see  below)  so  that  this  is  the  most  unstable 
inviscid  mode  for  a  hypersonic  boundary  layer.  In  Figure  (3.2)  we  show  the  eigenfunction 
of  the  vorticity  mode  equation  at  different  values  of  the  wavenumber.  This  figure  indicates 
that  as  the  wavenumber  decreases  the  eigenfunction  starts  to  expand  out  of  the  transition 
layer.  We  note  here  that  c<  as  defined  above  is  independent  of  Tj  and  that  the  growth  rate 
is  obtained  from  (3.3b)  by  dividing  dCj  by  +  1)(7  —  !)•  It  follows  that  wall  cooling 

has  a  destabilizing  effect  on  the  vorticity  mode. 

The  structure  of  the  growth-rate  curves  at  small  wavenumbers  is  of  interest  because 
for  sufliciently  small  values  of  the  wavenumber  we  expect  that  the  vorticity  mode  will 
develop  a  structure  similar  to  that  of  the  acoustic  mode.  In  the  following  discussion  we 
will  isolate  the  different  significant  regimes  which  occur  in  the  small  wavenumber  limit. 
For  simplicity  we  shall  now  concentrate  on  the  adiabatic  problem  and  take  Ti,  =  1. 

The  key  to  understanding  the  subsequent  regimes  when  a  is  related  to  inverse  powers 
of  the  Mach  number  is  to  write  down  the  small  a  asymptotic  structure  of  (3.4).  Figure 
(3.3)  is  a  schematic  illustration  of  the  different  regions  in  fj  space  wliich  emerge  in  the  limit 
a  — >  0.  Also  shown  in  this  figure  is  the  wall  layer  in  which  the  temperature  becomes  of 
order  M^.  For  the  moment  d  is  not  considered  to  be  sufficiently  small  for  d  to  be  0{M~*) 
for  some  positive  (j)\  it  then  turns  out  that  the  regions  I  and  IV  .shown  in  figure  (3.3)  are 
quite  passive.  However,  at  sufficiently  small  values  of  d  the  wall  layer  structure  of  the 
basic  state  will  enter  the  problem  -  see  below. 


After  some  careful  numerical  calculations  at  small  values  of  d  we  deduced  that  c 


expands  in  the  form 


c  = 


Cl  C2 

^4/7 


+  .... 


(3.9) 


Also,  from  (3. 2b, 4)  it  follows  that  for  ~  the  pressure  decays  like  exp(— d/;).  This 
suggests  that  the  pressure  in  III  should  be  expanded  in  the  form 


p  =  1  +  a^pi  +  d^p2  -h  ■  -h  a^P6  +  d(p7  -  T})  + 
d^(p8  —  Pit?)  +  d’ (p9  —  P2f))  +  d^(pio  —  P3t?)+  (3.10) 

+  d'^Pii  +  d^ Pi2  +  d^ Pi3  +  d^Pi4  +  .  . . , 


where  we  have  anticipated  the  form  of  several  terms  in  this  expansion  in  that  pi,  P2)  •  •  •  >  PlO 
are  taken  as  constants,  whilst  Pii,Pi2-,  etc.  are  functions  of  f).  At  order  dV"  in  this 
expansion  we  find  that  Pn  satisfies 


2G' 

Cl 


=  0, 


so  that  after  use  of  the  exponential  mathching  condition  it  follows  that 


2 

i’ll--—/  Gdj?-P4^  +  pii, 
Cl  U 


(3.11) 


where  pn  is  another  constant.  ^2  and  Pu  satisfy  similar  equations  with  forcing  functions 
coming  from  the  higher  order  terms  in  the  expansion  of  the  wavespeed. 

The  order  d^  term  in  (3.10)  is  then  found  to  satisfy 


/t  ^\2  2G  ,  C2.  ,C3  C2..  C4  2c2C3  cl. 

Pl4  =  (1  -  Gy  -  — (P3  -  — P2  -  (—  -  — )P1  -  —  +  -T2 - 

Cl  Cl  Cl  Cl  Cl  cf  Cl 


which  may  be  integrated  twice  to  give  P14.  However  it  is  enough  for  our  purposes  to  note 
that  when  f)  —*  0,  Pi4  ~  so  that  as  f?  — >  0  the  order  d^  and  d^  terms  in  (3.10)  become 
comparable  when  fj  =  0{ar)]  thus,  as  anticipated  earlier,  the  sublayer  III  is  of  depth  d^ . 
Within  this  layer  we  define 

Z  =  d"H 


and  expand  the  pressure  as 


p  =  1  +  d’pi  +  . . .  +  dp7  +  a^Ps{Z)  +  . . . , 


where  we  have,  by  matching  with  the  solution  in  II,  again  anticipated  several  terms  in  this 
expansion.  At  order  d’  we  find  that  Ps  satisfies 


fi//  72  81 

®  Z(ciZ4  +  9)^®  “ 


We  write  the  solution  for  Ps  in  the  form 


81dZ  _  1  . 
(ciZ'^  +  9)2 


where  the  constant  of  integration  has  been  chosen  to  satisfy  Pg  — v  —  1  as  Z  ^  oo  in  order 
to  match  with  the  solution  in  11.  Finally  we  consider  Region  IV  and  define  fj  =  a^z.  The 
zeroth  order  approximation  to  (3.4)  here  is 


8  81 

P»z  +  -Pz - gP  =  0, 


(3.12) 


which  must  be  solved  subject  to  p  — »  0  as  2  — ♦  0.  The  appropriate  solution  then  has  p 
constant  as  2  00  so  a  match  with  III  can  only  be  achieved  if 


p  81<iZ  1 

<0  +  ■ 


(3.13) 


It  is  assumed  that  ci  in  (3.13)  is  complex  so  that  after  a  little  manipulation  we  obtain 


which  corresponds  to  an  unstable  mode  of  (2.12a). 

g  =  0{M-^) 

Next  we  consider  the  situation  when  a  is  so  small  that  IV  in  Figure  (3.3)  merges  with 
the  wall  layer  of  the  basic  state.  At  this  stage  the  wall  layer  is  of  thickness  M~  >  so  that 
a  ~  M~  » .  The  wavenumber  a  is  then  written  in  the  form 

go  , 

g  = - r  +  . . . , 

Mf 

where  for  the  moment  we  shall  not  be  precise  about  the  size  of  the  first  correction  term  in 
this  expansion.  The  zeroth  order  approximation  to  (2.12a)  in  the  wall  layer  then  becomes 

^ ^  =  0.  (3.14) 
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where  a  dash  denotes  a  derivative  with  respect  to  the  wall  layer  variable  ^  and 

fZo  is  the  first  term  in  the  expansion  of  G  in  that  layer.  The  above  equation  is  again  to  be 
solved  subject  to  p'(0)  =  0.  For  large  ^  it  has  the  asymptotic  solutions 


p  ~  iVo  =  constant  and  p(  ~  iVi 


=  constant. 


For  most  values  of  oq  the  constant  Nq  is  nonzero,  and  the  structure  in  layers  I,  II,  III 
survives  intact.  Thus  for  these  values  of  ao  the  wavespeed  c  expands  as 


c  =  1  — 


(3.15) 


and  di  will  be  complex  so  that  the  wave  growth  rate  is  of  order  M~ 

However,  equation  (3.14)  has  a  countable  infinite  set  of  eigenvalues  for  which  the 
constant  Nq  =  0.  In  this  case  the  eigenfunction  takes  on  its  greatest  value  in  the  wall 
layer,  i.e.  where  the  steady  velocity  field  is  adjusting  to  it’s  free-stream  value.  These 
eigenvalues  correspond  to  the  acoustic  inviscid  modes  and  are  the  counterpart  of  those 
discussed  by  CH.  A  numerical  solution  of  (3.14)  yielded  the  sequence  ao  =  2.47,  7.17, 

12.19,  17.33,  22.54,  27.79,  _ The  first  three  eigenfunctions  associated  with  this  sequence 

of  eigenfunctions  are  shown  in  Figure  (3.4).  Thus  at  a  countable  discrete  set  of  points  the 
acoustic  modes  emerge  as  almost  neutral  “vorticity”  modes.  In  fact  we  could  seek  acoustic 
modes  with  a  ~  0(1),  i.e.  for  a  wavenumber  much  larger  than  currently  assumed,  for 
which  the  eigenfunction  is  concentrated  in  the  ^  =  0(1)  wall  layer.  These  modes  have  a 
fast  variation  in  this  layer  which  may  be  described  by  the  WKB  method.  At  certain  values 
of  M  these  eigenvalues  coalesce  with  the  neutral  vorticity  mode  discussed  earlier,  and  an 
analysis  outlined  in  CH  (see  also  SB)  can  be  performed  to  describe  the  ‘splitting’  of  the 
eigenvalues  in  this  region,  we  do  not  pursue  that  calculation  here  but  simply  determine 
the  structure  of  the  vorticity  mode  in  the  vicinity  of  the  acoustic  mode.  Thus  the  small  a 
structure  outlined  above  remains  valid  for  a  =  0{M~  > )  if  the  acoustic  modes  are  avoided; 
it  follows  that  continues  to  grow  in  the  further  limit  aM  >  — »  0,  where  no  neutral  acoustic 
modes  exist. 

In  order  to  examine  the  interaction  of  the  vorticity  mode  with  the  acoustic  modes  in 
a  little  more  detail,  we  expand  the  wavenumber  in  the  form 


jl/3/2  ^  25/14  ^ 
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whilst  c  still  expands  as  in  (3.15).  Here  we  are  denoting  by  oq  eigenvalue  of  (3.14) 
corresponding  to  No  =  0.  In  the  wall  layer  p  expands  as 

P=m2/'^Qo(H-^(M])  +  Qi  +  ...,  (3.16) 

where  we  have  anticipated  the  possible  occurrence  of  further  eigenfunctions  between  Qo,  Qi 
by  inserting  the  factor  (1  +  $[M]).  For  our  purposes  it  is  not  necessary  for  us  to  calculate 
$  here.  If  we  substitute  the  above  expansion  into  Rayleigh’s  equation  we  see  that  Qq 
satisfies  (3.14)  with  a  =  Oq.  The  derivative  of  Qq  vanishes  at  ^  =  0  whilst  for  ^  >>  1 

Qq  ~  qoC^ 


where  qo  is  a  constant.  At  higher  order  in  Rayleigh’s  equation  we  find  that  Qi  satisfies 

„  201,  V  (ao)^(^-l)^,  .2,,,  2(1-10) 

Ql  -  — - tQi - : - -  “o)(l  -  77 - 7T7TT7r-T)vl 


Go  —  1 


4  '  (^_1)(1  +  Q0) 

- 1)2(1  -  a§)(l  - 


2  '•  '  (7-l)(l  +  ilo) 

The  solution  of  this  equation  which  has  Qi(0)  =  0  is  such  that  for  large  ( 


Ql  ~  ooqooi* 


The  constant  ao  depends  on  oq  and  7,  and  can  only  be  determined  numerically.  However  it 
can  be  shown  that  uq  is  alternatively  positive  and  negative  at  successive  values  of  Oq.  The 
region  III  of  Figure  (3.3)  becomes  of  depth  0(M“^/^^)  when  a  is  The  solution 

in  this  layer  is  calculated  using  the  procedure  outlined  above;  the  only  significant  change 
is  that  in  matching  regions  III  and  IV  we  must  now  account  for  the  fact  that  Qo  ~ 
for  large  After  some  manipulation  we  find  that  the  eigenrelation  obtained  when  this 
matching  has  been  carried  out  is 


J_  -  J_  1 
“1  ~  if  ^ 


(3,17) 


where 


«l  = 


BA^aodo^l 

7D^ 


d'l  =  Ddi,  D  =  (9— A  = 
Sy/2 
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(7  -  1) 


(3.18a,6,  c,  d) 


In  Figure  (3.5)  we  have  shown  solutions  of  (3.17)  for  — 5  <  <  5.  At  both  ends  of  this 

Figure  di  approaches  e~  ^  which  gives  the  required  matching  with  the  unstable  vorticity 
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modes  away  from  the  acoustic  modes.  We  further  note  that  for  small  the  growth  rate 
tends  to  zero.  More  precisely  from  (3.17)  we  find  that  for  small  aj 


a 


%l/8 

\j*  A 


r/4 


+ 


SO  that,  dependent  on  on  the  sign  of  aj,  the  growth  rate  goes  to  zero  like  M~  >  or 

Strictly  we  should  now  rescale  the  above  expansions  when  aj  becomes  sufficiently 
small.  We  do  not  pursue  that  rescaled  problem  here  since  our  main  concern  is  with  com¬ 
pleting  a  discussion  of  the  structure  of  the  vorticity  mode  at  all  lengthscales,  in  particular 
we  wish  to  see  if  the  vorticity  mode  connects  with  any  other  neutral  states.  We  have 
already  seen  that  Cj  continues  increasing  as  the  wavenumber  decreases  through 
values,  apart  from  small  neighborhoods  of  the  acoustic  modes  where  it  decrease  towards 
zero.  Thus  Cj  tends  to  infinity  on  the  present  scaling  when  ao  tends  to  zero.  Thus  the 
main  significance  of  the  a  ~  regime  is  that  it  is  at  this  stage  where  the  acoustic 

modes  emerge;  however  the  small  a  structure,  developed  initially  for  the  a  =  0(1)  scaling, 
survives  this  regime  intact  in  the  limit  of  ao  tending  to  zero. 


g  =  OjM-i) 

The  next  significant  stage  in  the  development  of  the  vorticity  mode  occurs  when  the 
temperature  in  the  upper  layer,  i.e.  where  ~  g~^,  becomes  such  that  T  ~  (tz  —  c)^M^. 
The  pressure  eigenfunction  in  the  upper  layer  then  has  its  decay  to  zero  modified.  This 
situation  occurs  when 


with 


(3.18) 


Q  = 


M'^/4 


+ 


(3.19) 


In  this  case  the  pressure  in  the  upper  layer  decays  like  exp(— f(l  —  df)*  This 

means  that  in  the  rescaled  region  corresponding  to  III  of  Figure  (3.2),  where  we  write 
Z  =  the  pressure  perturbation  expands  as 


p  =  1  + 


P2{Z) 

m2 


and  P2  satisfies  the  equation 


p"  _ 

^  Z[diZ^-3A]  ^  Z8 
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The  solution  of  this  equation  which  enables  us  to  match  with  the  behaviour  discussed 
above  is  _ 


f.,  ,1  t- 

P2  -  'V 

In  the  wall  layer  there  are  now  no  acoustic  mode  eigenfunctions  which  lead  to  decay  at 
infinity  so  that  matching  leads  to  the  eigenrelation 


9^2(7  -  1)2 
(diZ^  -  3A)2 


dZ  + 


l\/ 1  —  (if 


A  more  convenient  form  for  the  eigenrelation  is  found  by  writing 


8ls/2 

3(3A)^/‘*(7  —  l)27r 


so  that 

As  expected  in  the  limit  £  —*  oo  we  obtain  di  ^  g-3<ir/7£-  i 

which  is  the  limiting  form  of 

the  solution  obtained  for  a  ~  M“  > .  It  is  a  simple  matter  to  show  from  the  above  equation 
that  du  is  a  monotonic  decreasing  function  of  £  and  in  particular  for  f  »  0  we  find  that 


(3.20) 

The  main  implication  of  this  limit  is  that  the  M~  *  layer  for  the  a  ^  M~^  modes  decreases 
in  size  like  *  when  £\  —*  0.  The  next  stage  is  when  the  M~*  and  M"  »  layers  merge 
which  occurs  when  di  =  0{M)  and  a  = 

a  =  O(M-t) 

Here  we  expand  c  and  a  in  the  forms 

L 

c  =  CQ  +  . . . ,  a  =  — T-  + - 

Mi 

Now  there  are  just  three  regions,  of  depth  M~  >  ,M®,  and  M» ,  for  us  to  consider.  Inter¬ 
estingly  the  extent  of  the  perturbation  away  from  the  wall  has  been  reduced  by  a  factor 
of  M«  from  the  order  M«  scale  of  the  case  a  ~  M"  > .  This  reduction  in  depth  is  caused 
by  (c  —  c)  remaining  0(1)  in  this  layer,  and  if  we  write  rj  =  M^Y,  the  dominant  terms  in 
the  equation  for  the  pressure  in  that  layer  give 


Pyy  +  -  co)^p  =  0 


(3.21) 
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If  CO  has  a  positive  imaginary  part,  it  follows  that 


p  =  exp(-(Ty), 


with  <7  =  iL{l  —  cq).  Thus  where  y  — »  0  we  have  that  — »  —M  ’<t.  This  means  that 

the  required  solution  for  p  in  the  t]  =  0(1)  region  is 


p  -  1  +  +  ^]  + -  (3.22) 

Here  the  constant  terms  are  small  compared  to  but  large  compared  to  M~  in  prin¬ 
ciple  they  may  be  found  by  matching  requirements.  It  remains  to  consider  the  small  layer 
when  T)  =  M~  3^,  and  p  expands  as 


P='^+---  +  -^^2iO  + 


p2  is  found  to  satisfy 


^2  -  -  1)(1  +  wo)(l  -  wo){^(7  -  1)(1  +  wo)(l  -  uo)  -  (uo  -  cf) 

uo  —  CO  2  2 

SO  that 


o/  _  ,2  '^^(7-l)(l  +  wo)(l-Uo),(7-  l)(l  +  uo)(l-uo) 

2^2=(“0-Co)  /  - - 72 - ( - r - 

Jq  2(uo-co)^  2 

A  match  with  the  core  solution  is  achieved  if 


(uo  -  co)^)d^ 
(3.22) 


L2(7_i)(i  +  uo)(1-uo),,  ,,,,  ,  ,,,  , 

- -  co)^ -  ^  ^  ^  °  ^  ~  ~  ^  °  ~ 

=  -cr  =  -iL(l  -  co). 


(3.23) 

This  is  the  required  eigenrelation  to  determine  the  complex  wavespeed  co  as  a  function 
of  the  scaled  wavenumber  L.  A  large  L  analysis  of  this  equation  shows  that  the  limiting 
small  £  solution  of  the  a  =  0(M~  *  )  case  is  retrieved. 

In  Figure  (3.6)  we  show  co  as  a  function  of  L  for  100  <  L  <  400  and  the  case  of  an 
insulated  wall.  We  see  that  the  maximum  value  of  c^  for  a  two-dimensional  mode  occurs 
for  L  ~  225  and  a  neutral  mode  exists  for  L  ~  110.  We  postpone  until  the  final  section  a 
discussion  of  the  implications  of  the  results  we  have  found  above. 
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4  The  Inviscid  Instability  Problem  in  Interactive  Boundary  Layers. 

We  consider  the  hypersonic  flow  of  a  Sutherland  law  fluid  past  the  semi  infinite  wall 
defined  by  the  positive  x  axis.  The  oncoming  flow  has  constant  velocity  at  infinity  and 
the  leading  edge  of  the  wall  is  sharp.  This  results  in  the  formation  of  a  shock  wave  which 
acts  as  an  upper  bound  for  the  flow  disturbances  in  the  sub-layers  below,  the  uppermost 
of  which  is  inviscid.  The  reader  is  referred  to  Figure  (4.1)  for  a  description  of  the  different 
parts  of  the  flow  field. 

This  region  has  been  studied  by  for  example,  Bush  (1966),  and  since  viscosity  is 
negligible  here  the  use  of  a  power- viscosity-temperature  law  does  not  alter  the  well-known 
governing  equations  for  this  region.  Below  the  expansions  and  equations  governing  the 
flow  are  given.  Apart  from  some  minor  differences,  our  formulation  is  essentially  the  same 
as  that  of,  amongst  others,  Stewartson  (1955,  1964)  and  Bush  (1966)  and  so  the  reader  is 
referred  to  those  papers  for  more  details. 

The  flow  quantities  and  coordinates  have  been  non-dimensionalised  using  upstream 
flow  quantities  and  L  a  typical  streamwise  lengthscale.  We  then  scale  the  non-dimension¬ 
alised  normal  coordinate  y  appropriately 


y  = 
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M’ 


(4.1) 


where  M  is  the  Mach  Number  and  Y  is  taken  to  be  0(1).  Note  that  as  yet  we  do  not 
scale  the  downstream  variable  x.  As  in  previous  studies  we  consider  the  two-dimensional 
problem  and  this  enables  the  introduction  of  a  streamfunction  ip  defined  by 


PU  =  IPy  ,  PV  =  -1P^ 


(4.2a, 6) 


Further  we  expand  the  velocity,  temperature,  density,  pressure  and  streamfunction  as 


ui 

“=l+li72  + 


r  =  Ti  -t- 


P  =  PI  + 


PI  , 


(4.3a,  6,  c,  d,  e,  /  ) 


and 


^  =  77^1  + 
M 
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If  the  above  expansions  are  substituted  into  the  inviscid  Navier-Stokes  equations  written 
in  Von-Mises  coordinates  we  obtain,  after  some  manipulation,  the  governing  equations  for 
this  layer 

vu  =  -PU,  ,  vu,  ^  E{rpi)pl.  (4a, 6, c) 

The  function  Eitpi)  is  evaluated  from  the  initial  conditions  which  for  these  hyperbolic 
equations  are  given  at  the  shock.  Conventionally,  we  define  the  shock  by  y  =  f(x)  where 
/  is  unknown  at  present;  we  note  that  =  /  on  the  shock. 


The  Strong  Interaction  Zone 

The  solution  of  these  equations  (and  the  corresponding  ones  for  the  lower  layers)  can 
be  investigated  analytically  in  the  large  x—  and  small  a:— limits  using  expansion  procedures, 
following  Stewartson  (1955,  1964),  Brown  and  Stewartson  (1975).  The  formulation  below 
follows  the  latter  paper  closely. 

Let  us  consider  the  solution  of  (4.4a,b,c)  for  small  x,  i.e.  close  to  the  leading  edge  of 
the  plate,  where  the  shock  is  attached.  Here  we  write 

/  =  +  •  •  •  ,  Pi  =  ai?i(C)a;"  1/2  ^  ^ 

ui  =  aivi(C)x"^/‘^  +  •••  and  Pi  =  Pi(C)  +  •  •  •  .  (4.5a,  6,  c,d) 


The  similarity  variable  C  introduced  in  the  above  equations  is  defined  by 


aix^/4 


(4.6) 


The  scalings  (4.5b, c,d)  are  implied  by  the  Rankine-Hugoniot  relations,  which  relate 
flow  quantities  either  side  of  the  shock  Y  =  f{x),  and  by  (4.5a)  which  follows  from  matching 
normal  velocities  and  pressure  across  the  sub-layers  beneath  the  shock;  see  Stewartson 
(1964)  and  references  therein. 

The  resulting  equations  for  the  leading  order  terms  are 

+  ^ClTi  =  4^1, 

^  9(7-1)'' 

=  and  ei  =  (4.7a,  6,  c,d) 

with 

3  _  9  _7+l  ,, 

"‘=2(TTTj’  ’’*  =  5(7TTj’  ^  = 


20 


which  are  the  conditions  at  the  shock.  We  require  their  solution  for  C  <  1-  In  particular, 
for  given  7  we  can  solve  numerically  for  ui(0)  and  Pi(0). 

Substitution  of  (4.5b,d)  into  (4c)  and  recalling  the  definition  of  (  gives  E{'tpi)  ~ 
PlOi  xpi  as  V’l  0.  This  implies  that,  as  ^  >  0 ,  ~  and  this  produces  a 

higher  order  correction  in  ui  if  7  >  For  a  realistic  gas  7  is  larger  than  this  value, 
with  the  result  that  there  must  be  a  viscous  sublayer  beneath  the  present  layer.  Bush 
(1966)  found  that  this  layer  enables  one  to  match  the  inviscid-layer  solution  on  to  the 
boundary  layer  solution  if  one  used  a  power-viscosity/ temperature  law  (p  ^  T“  ,u>  <  1), 
rather  than  the  linear  law  used  by  Stewartson  and  most  other  research  studies  in  the  field. 
Stewartson  (1955)  hoped  that  by  using  idealized  assumed  physical  properties  it  would 
help  to  “understand  the  behaviour  of  more  realistic  fluids”.  His  assumption  resulted  in 
discontinuous  derivatives,  or  “kinks”,  in  properties  between  the  Inviscid  Zone  and  the 
Boundary  Layer,  see  Fig.  7.4,  page  167  of  Stewartson  (1964).  Despite  the  results  of  Bush 
(1966)  (and  a  few  subsequent  authors)  the  linear  viscosity  law  is  still  extensively  used 
in  hypersonic  shock/boundary-layer  research;  it  is  argued  that  these  ‘kinks’  are  only  a 
slight  nuisance  that  can  probably  be  explained  away  by  deeper  analysis  and  appeal  to  the 
argument  of  Stewartson  concerning  the  assumption  of  ideal  fluid  properties.  The  analysis 
of  Bush  was  extended  to  the  linear  case  by  Lee  and  Cheng  (1969).  They  found  that,  by 
consideration  of  a  second  order  boundary-layer  correction,  coupled  to  Bush’s  transitional 
layer  analysis,  the  two  layer  structure  of  Stewartson  gives  the  correct  basis  for  the  flow 
structure. 

Not  surprisingly  the  extensive  use  of  Chapman’s  law  in  hypersonic  calculations  has 
encouraged  researchers  in  stability  theory  to  make  the  same  rather  severe  approximation  to 
the  fluid  viscosity;  we  believe  that  the  calculations  given  in  this  paper  are  the  first  to  take 
account  of  a  realistic  viscosity-temperature  dependence  in  a  first  mode  stability  analysis. 
We  note  that  our  flow  solutions  in  the  layers  beneath  the  shock  are  in  full  agreement  with 
Bush,  but  he  considered  a  general  power  law  whereas  we  concentrate  on  p  the 

leading  order  form  for  Sutherland’s  law  at  high  temperatures.  For  convenience  we  will 
use  the  notation  of  Stewartson  and  Brown  (1975).  We  will  next  consider  the  viscous- 
adjustment-layer  after  discussing  the  third  sub-layer,  the  viscous  boundary  layer. 

Viscous  Boundary  Layer 

This  is  adjacent  to  the  surface  of  the  plate  and  enables  velocities  to  be  reduced  to 
zero  on  the  surface.  First,  the  full  asymptotic  formulation  is  described  briefly,  then  the 
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asymptotic  behaviour  for  small  x  is  discussed  in  the  next  subsection.  We  assume  that  the 
boundary  layer  exists  where  V’  ~  M~  ,  V’  being  the  stream-function  defined  earlier.  We 
define  the  hypersonic  parameter  r  by 


Re  =  rM^ 


and  take  r  to  be  order  one.  Note  that  the  power  of  the  Mach  number  is  six  in  the  definitions 
of  the  hypersonic  parameters  of  those  using  a  linear  viscosity  law.  We  write 


u  =  Ui-{- 


^  =  77  + 
M 


T  =  M^6\  and  n  =  Mm. 


{4.9a,b,c,d,e,f) 


These  scalings  are  standard  apart  from  that  for  the  viscosity  n  which  follows  immediately 
from  the  high  temperature  form  of  Sutherland’s  law: 


M  =  (1  + 


(4.10) 


We  make  the  assumptions  that  the  Prandtl  number  is  unity  and  that  the  wall  is 
insulating,  this  enables  us  to  carry  out  a  relatively  simple  analysis  for  the  viscous  solutions. 
The  solutions  capture  the  new  structure  of  the  inviscid  modes  due  to  the  effect  of  using 
a  power-viscosity  law  in  the  viscous  sublayers.  The  relaxation  of  these  assumptions  does 
not  substantially  alter  the  conclusions  obtained  below. 

We  again  transfer  to  von  Mises  coordinates;  using  the  streamfunction  ^  as  an  inde¬ 
pendent  variable  instead  of  Y.  The  energy  and  momentum  equations  are  combined  to 
yield  the  boundary-layer  equations  in  these  coordinates 

2  7P1  r  \9i  y  ^ 


h  =  ^(1  -  t'l). 


1/2 

together  with  m  =  from  (4.10),  and  with  usual  boundary  conditions 


(4.11a,6) 


Ui  =  0  on  =  0, 


Ui  —y  1  as  ’J'  -+  00. 


(4.12a, 6) 
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Boundary-Layer  Solution  for  the  Strong  Interaction  Zone 

Now  let  X  be  small  so  that  the  hypersonic  interaction  parameter  becomes  large.  We 
define  the  similarity  variable  (f)  by 

/  o  \  t.1/4 

*  =  [47(1  +  C)p,{0)|l'2(^— )  (4.13) 

where  ai  and  Pi(0)  were  introduced  when  considering  the  inviscid  region.  If  we  now  write 

Ui=gi<i>)  +  ---,  (4.14) 


we  obtain  the  similarity  equation  for  the  boundary-layer  solution 

-to'  =  (^)  (1  -  s")  +  9(75-1-^) 


(1  -  ^2)1/; 


and 


i?(0)  =  0  ,  g{cx)  =  1  . 


(4.15a,  6,  c) 


This  is  a  modified  Faulkner-Skan  equation,  but  the  decay  of  g  at  large  4>  is  now 
algebraic  rather  than  the  usual  exponential  behaviour. 


a 


+ 


as 


00, 


(4.16) 


and  a  >  0  can  be  determined  analytically. 

Matching  normal  velocities  across  the  sub-layers  in  the  usual  manner  yields  the  leading 
order  coefficient  ai  in  the  small  x  expansion  for  the  as  yet  unknown  shock  location 


a? 


4ui(0)  N  r7Pi(0)  '  ^  2  ^  J Q  g 


(4.17) 


The  Viscous  Transition  Layer 

Here  the  velocities  and  temperatures  in  the  outer  inviscid  flow  match  with  those  in 
the  inner  boundary-layer  flow.  We  give  a  brief  account  of  the  arguments  used  to  deduce 
the  position  and  properties  of  this  region  and  again  consider  the  strong  interaction  zone. 

Earlier  it  was  noted  that  the  inviscid  density  function  ~  C  0-  It  can 

easily  be  deduced  that  ui  and  Ti  grow  like  inverse  oowers  of  C  as  C  — »  0.  In  fact 

/  ™3/4\  2/37 

r~ri  )  1  (4.18) 

\  ^1  / 
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as  V'l  — »  0. 

Now  let  us  consider  the  boundary-layer  lemperature  as  the  top  of  the  layer  is  ap¬ 
proached.  We  have  written  T  =  and  found  that  9i  =  -*^(1  —  U^)-  We  know  the 

large  (p  behaviour  of  17i  and  it  follows  immediately  that 

SO  that 

T  =  ~  —f  -  „■  as  p  oo.  (4.19) 

Matching  the  limiting  temperatures  (4.18)  (4.19)  gives  the  position  of  the  viscous  transition 
layer, 

,  A  =  -  (4.20a,  6) 

67  —  1 

This  is  in  full  agreement  with  Bush’s  result.  Note  that  the  powers  involved  in  the  asymp¬ 
totics  are  functions  of  7,  the  ratio  of  specific  heat  capacities,  in  the  new  layer.  Further 
investigation  of  the  small  pi  limit  of  the  inviscid  solutions  implies  the  following  expansions 
for  flow  properties  in  the  transition  layer 

u  =  1  +  M  ui  4-  •  •  •  , 

r  =  MJ^ri  +  ---,  (4.21a,6,c) 


together  with 


P  = 


7P1 

T  ' 


We  define  the  scaled  streamfunction  77 


=  -qM  *  (=  Mp). 


(4.22) 


These  are  then  substituted  into  the  Parabolised  Navier-Stokes  equations  in  Von-  Mises 
co-ordinates  to  obtain  the  equations  for  ui 

uu  ^  Tipi,  iPid  (  1  auiA 
(1  -f  C)  -fpx  r  dq\  jk  dq  r  ^  > 

where,  as  we  are  taking  the  Prandtl  number  to  be  unity  and  the  wall  to  be  insulating,  the 
corresponding  energy  equation  can  be  integrated  to  yield 


Ti  =  (1  -  7)«1- 


(4.24) 
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The  pressure  pi  is  independent  of  the  normal  co-ordinates,  as  usual  in  viscous  regions 
governed  by  the  boundary-layer  equations,  so  the  pressure  is  expanded 

fllPl(O)  , 

Pi  =  - 1 —  +  •  •  • 

a; » 

in  the  two  viscous  layers,  this  being  the  limiting  form  as  the  inviscid  layer  is  descended. 

We  see  that  (from  (4.23))  in  this  adjustment  layer  we  still  have  the  high  temperature 
form  of  Sutherland’s  law,  in  contrast  with  the  shock-free  adjustment  layer  of  section  2 
where  the  full  form  is  needed.  Again,  a  similarity  solution  for  the  flow  exists  in  this 
viscous  sub-layer  with  similarity  variable  z  defined  by 

T)  =  /3x^z,  (4.25a) 


where 

^2  ^  (C+  l)7ail>l(0) 
r 


and 


Writing  ui 
crucial  region 


A  = 


3-A 


(4.256, c) 


x^~^G(z)  leads  to  the  governing  equation  for  the  flow  solution  in  this 


(A  -  2)G  -  ^ — ^zG'  =  — -G  + 
4  27 


V{l-l)G 


(4.26) 


Note  that  from  the  definition  of  ui,  (4.21a),  we  expect  that  G  <  0.  We  require  that  this 
solution  for  ui  matches  to  the  corresponding  solutions  above  and  below,  i.e.  in  the  inviscid 
sub-layer  and  viscous  boundary  layer  respectively,  so  we  must  investigate  asymptotically 
the  small  and  large  z  limiting  forms  of  G. 

As  ir  — »  0  we  find  that  G  has  the  behaviour 


G  =  -goz~'^  +  giz  ’'  +  •••, 


_ 5767^ 

(7  -  1)(37  -  1)2’ 


(4.27a,  6) 


Thus  the  transition  layer  solution  matches  onto  the  large  (j)  form  of  g.  The  exponent  u 
of  z  in  the  first  correction  term  can  be  calculated  analytically  whilst  the  coefficient  gi 
must  be  determined  from  a  numerical  solution.  In  fact  u  satisfies  a  quadratic  equation 
with  coefficients  functions  of  7.  With  7  =  1.4  we  find  that  1/  =  0.6267 ...  so  that  we  have 
relatively  small  correction  terms.  The  value  of  go  for  this  choice  of  gamma  is  275.6  . . .  and 
so  G  grows  quickly  as  i  — >  0. 
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We  now  investigate  the  behaviour  of  G  for  large  z  we  find  that 


G  —*■  —Aqz  ^  +  Aiz  ^  +  •  •  • , 


(4.28) 


and  so  again  we  have  a  relatively  small  correction  term.  Note  that  G  decays  algebraically,  in 
contrast  to  the  rapid  exponential  decay  of  the  Blasius  and  the  “Modified  Blasius”  functions 
which  arise  in  the  shock-free,  far-downstream  cases  (the  former  after  employing  the  linear 
Chapman  viscosity  law  and  the  latter  from  the  use  of  the  more  realistic  Sutherland’s  law 
-  see  earlier).  The  second  coefficient  Ai  is  a  function  of  gamma  and  can  be  determined  in 
terms  of  Aq;  in  fact 

84(37  +  1) 

97(7  —  !)>  (37  —  1) 

The  constant  Aq  is  not  fixed  by  the  asymptotics.  Instead  we  must  choose  it  so  that  the 
velocity  solution  matches  with  the  appropriate  limiting  form  of  the  inviscid  solution.  We 
obtain,  after  an  elementary  matching  argument,  and  a  little  manipulation,  the  result 


Ao  =  (1  +  C)- 


_l_  o  J_  1- 


4 

*7 


1 


7  -  1 


(4.29) 


For  the  value  7  =  1.4  we  find  Ao  =  0.1751....  In  our  numerical  calculations  we  took 
without  loss  of  generality  the  hypersonic  interaction  parameter  to  be  unity.  Now  that  the 
base  flow  for  this  region  has  been  found  (at  leading  order)  we  can  consider  its  stability 
characteristics;  in  particular  we  are  interested  here  in  the  linear  stability  of  inviscid  modes 
concentrated  (trapped)  within  this  adjustment  layer.  These  so-called  vorticity-modes  have 
been  discussed  earlier  for  other  flows. 


The  Vorticity  Mode  in  the  Strong  Interaction  Zone 

The  scalings  for  these  modes  appear  complicated  but  follow  in  a  straight-forward 
manner  after  applying  the  usual  vorticity  mode  arguments  to  the  flow  field  discussed 
above.  The  modes  of  wavenumber  k  have  wavelength  comparable  with  the  thickness  of 
the  transition  region  and  a  Rayleigh  analysis  suggests  we  require  9^^  ~  This  can 

be  seen,  for  example,  from  (2.12a)  and  (3.4)  studied  earlier.  From  (4.21b),  and  recalling 
the  definition  of  A  in  (4.20b),  we  have  T  and  also  from  (4.22),  xj)  M 

Note  that  with  these  scaling  we  are  primarily  concerned  with  obtaining  the  correct  Mach 
number  dependence;  the  small  x  dependence  will  be  incorporated  later.  We  deduce  that 


26 


k  ^  >>  1,  and  thus  this  represents  a  short  wavelength  mode  and  so  we  introduce 

a  fast  x-scale 

X  =  (4.30) 

The  time-scale  can  also  be  similarly  deduced  from  Rayleigh  analysis,  after  noting  that  the 
vorticity  modes  propagate  in  a  frame  moving  with  uniform  velocity  and  recalling  the  ui 
scaling  (4.21a)  .  We  are  lead  to  the  introduction  of  a  short  time  scale. 

T  =  M^-h.  (4.31) 

We  note  that  these  scales  are  compatible  with  our  earlier  analysis  of  the  non-interactive 
case.  At  leading  order  we  have  the  multiple-scales 

d,  -^d^+  ,  d,  M^- -  M'^-^^dx  . 


Note  that  non-parallel  effects  are  0(3,,)  ~  0(1)  and  are  negligible  in  comparison  with  the 
direct  growth  effects  of  order  k(u  —  1)  =  >>  1;  for  our  choice  of  7  =  1.4  we 

find  that  A  =  1.7027.  The  wavespeed  c  of  the  vorticity  mode  propagating  in  the  frame 
moving  with  uniform  (nondimensionalised)  velocity  u  =  1  is  related  to  its  stationary-frame 
wavespeed  c  by 


c  =  1  + 


c 

A/4(2-»)  +  ■■■ 


(4.32) 


which  follows  immediately  from  (4.21a).  Note  that  it  represents  a  small  perturbation  about 
the  velocity  of  the  propagating  frame.  The  remainder  of  the  analysis  follows  the  classi¬ 
cal  inviscid  mode  approach  for  formulating  the  pressure  equation  for  the  linear  wave-like 
disturbance,  now  that  the  scales  have  been  deduced.  We  assume  that  the  flow  previously 
discussed  for  this  region  is  perturbed  by  an  infinitesimal  disturbance  proportional  to 


^ii(X-eT) 


(4.33) 


containing  all  the  X  and  r  dependence  of  the  disturbance.  We  note  that  the  basic  flow 
properties  are  independent  of  r  and  X  allowing  such  an  assumption.  After  a  little  manip¬ 
ulation  we  will  obtain  the  so-called  vorticity-made- Rayleigh  equation  for  the  amplitude  of 
the  disturbance  pressure  p,, 


PtlT) 


2ui,  ^  (7-1)M,.2^  n 

■P'1 - 7^ - ^  P  =  0, 


ui  —  c 


l^Pl 


(4.34) 
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where  we  require  p  to  decay  to  zero  for  large  and  small  t).  We  now  return  to  the  Strong 
Interaction  Zone  and  investigate  the  asymptotic  behaviour  of  (4.34)  for  small  x. 

We  recall  the  leading  order  small  x  forms 


afpi(O) 


PI  = 


ui  =  X 


X  a 

x-2 


TJ  =  0X*^^ 


G. 


(4.35a,  b,  c) 


We  need  the  small-x  dependence  of  c  and  k  but  need  not  worry  about  p  as  (4.34)  is  linear. 
We  find  that 


c  =  x^-^C, 


and 

=  (^_l)j3(A-l)/4^’  (4.36a,  6) 

are  the  necessary  behaviours.  The  former  follows  immediately  from  (4.35c)  and  the  latter 
results  from  considering  the  coefficient  of  pin  (4.34).  We  finally  arrive  at  the  renormalised 
problem  for  p;  the  vorticity-mode  pressure-amplitude  equation  for  the  strong-interaction 
zone 

Pzi  -  -  IC^G^p  =  0,  (4.37) 

to  be  solved  subject  to  p  vanishing  in  the  limits  z  0  and  z  — »  oo,  so  that,  as  in  previous 
analyses,  the  disturbance  is  confined  to  the  adjustment  layer.  We  discuss  our  numerical 
solution  of  this  eigenvalue  problem  and  its  limiting  small  wavenumber  form  in  the  next 
section. 
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5  The  Solution  of  the  Strong  Interaction  Vorticity  Mode  Equation. 
The  leading  order  asymptotes  of  equation  (4.37)  are  found  to  be 


PoZ  exp  I 


(  n 

,  .  _i-  (  37A0X:  (37-2)\ 


and  p  — »  Poo  exp  l  — ^ - -z- — r - -I  a,sz 00 

(5.1a, b) 

where  po  and  Poo  are  constants  and  the  numbers  50  and  Aq  arise  from  the  corresponding 
asymptotes  for  G  given  earlier.  Higher  order  terms  in  these  expressions  can  be  found 
analytically,  and,  in  fact,  are  needed  for  accurate  numerical  solutions  to  this  eigenvalue 
problem  for  C  and  /C,  the  scaled  relative  wavespeed  and  the  scaled  wavenumber  of  the 
inviscid  vorticity  mode,  discussed  after  briefly  paying  attention  to  the  neutral  case.  There 
C  is  real  and  takes  the  value  of  G  evaluated  at  the  generalized  inflection  point  which  of 
course  occurs  where 

=  2GI 


A  numerical  solution  to  (4.26)  using  a  Runge-Kutta  method  shows  that  this  occurs 
when  ^  =  1.661432,  where  the  new  variable  ^  =  in z  was  introduced  to  ‘stretch’  the  co¬ 
ordinate-scale  for  small  z  where  G,  and  hence  p,  vary  rapidly.  The  resulting  neutral  value 
of  the  wavespeed  is  then  given  by 


C  =  -0.633318. 

The  corresponding  real  value  of  K  is  obtained  from  a  numerical  solution  of  (4.37)  for  all  z 
with  an  appropriate  treatment  of  the  path  of  integration  at  the  generalised  inflection  point 
(g.i.p.).  The  particular  numerical  procedure  chosen  to  calculate  this  neutral  wavenumber, 
and  the  complex  eigenvalues  C  was  to  extend  the  line  of  integration  into  the  complex- 
plane,  taking  a  triangular  contour  around  and  below  the  g.i.p.  The  neutral  value  of  K 
was  calculated  to  be 

K  =  .477957 

Figure  (5.1)  shows  the  growth  rate  KCi  plotted  against  K,  we  find  that  the  maximum 
growth  rate  ~  .060918  occurs  at  /C  ~  .156100  and  that  the  growth  rate  goes  to  zero  when 
the  wavenumber  goes  to  zero.  In  Figure  (5.2)  we  show  the  eigenfunction  of  the  most 
unstable  mode,  we  note  the  exponential  decay  of  the  eigenfunction  at  both  ends  of  the 
r2uige  of  integration. 
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We  shall  now  show  how  the  small  wavenumber  structure  of  the  strong  interaction 
inviscid  mode  develops  when  the  wavenumber  tends  to  zero.  In  fact  we  simply  give  the 
structure  corresponding  to  the  first  small  wavenumber  regime  discussed  in  Section  3.  Some 
discussion  of  the  regimes  encountered  at  even  smaller  wavenumbers  will  be  given  in  the 
next  section. 

The  small-/C  behaviour 

Numerical  solutions  of  the  pressure-amplitude  equation,  (4.37),  indicate  that  C  grows 
as  /C  — »  0;  exactly  the  same  behaviour  was  found  for  the  no  shock  problem  studied  in 
Section  3.  However,  in  the  latter,  the  leading  order  dependence  on  the  wavenumber  was 
indicated  from  the  numerical  solutions,  whereas,  in  the  present  case,  such  guidance  was 
not  possible  due  to  a  complicated  dependence  ofC  on  /C.  Instead,  the  equation  was  studied 
analytically  for  small  /C,  and  initial  investigations  suggested  that 


(fc2)-4e  _  !■ 

4e 


c  ~ 


as  A:  — ►  0.  Here  the  parameter  e  is  a  function  of  7  and  is  defined  by 


4e  =  1  -  — 
37 


(5.2) 


and  we  have  replaced  1C,  C  hy  k,  c  respectively  here  and  henceforth  in  this  sub-section. 

4 

It  follows  from  (5.2)  that  special  attention  is  required  for  values  of  7  close  to  7  =  - 

o 

and,  as  this  coincides  with  the  range  of  physically  relevant  values  of  the  ratio  of  specific 

heat  capacities,  we  solely  consider  the  case  of  lc|  <<  1  in  this  section.  We  note  that  the 

4 

suggested  form  of  c  is  defined  and  continuous  at  7  =  -  by  appealing  to  “L’Hopital’s  Rule”. 

O 

Below  we  outline  the  asymptotic  regions  and  corresponding  solutions  required  to  resolve 
the  leading  order  behaviour  of  c  as  A:  — »  0. 

When  z  is  very  large,  we  have  the  WKB  asymptote 


P~ 


exp  I 


37^0 

37-2 


kz 


and  we  note  that  for  smaller  z  this  form  breaks  down  in  a  ‘turning-point’  region  corre¬ 
sponding  to 

Z  k~  >>  1. 


In  this  region  we  define 

yz=k^z,  P  =  Po{y)  +  --- 


(5.3a,  6) 
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and,  using  the  large-2  asymptotes  for  G,  we  find  that  Pq  satisfies 

Poyy  -  Aly~'^‘‘pQ  =  0,  a  = 

J7 


(5.4a, 6) 


This  equation  for  Pq  has  an  analytic  solution  involving  the  modified  Bessel  Function 
K,  the  required  solution  being 


where 


Po  =  Dor]'"K^{-r]) 


T]  =  \AQy.y''  ^  V  =  \  —  2a  and  y.  =  — . 

2v 


(5.5a) 


(5.5b,  c,  d) 


Here  Dq  is  an  arbitrary  constant.  We  can  investigate  the  small  y  behaviour  of  Pq  by 
considering  the  limiting  forms,  for  small  arguments,  of  the  modified  Bessel  functions.  We 
find  that 

p.  =  Ml  +  +  -I 

sm{yir)  (I  -  y)  T(1  -  y)(l  -  y) 

as  y  — »  0,  the  ordering  of  the  second  and  third  terms  being  dependent  on  the  value  of  7; 
we  note  that  the  corresponding  powers  of  y  (the  indicial  roots  of  the  associated  Frobenius 
solution)  become  equal  when  7  =  -.  Moreover,  the  coefficients  of  these  terms  are  also 

4 

singular  when  7  has  this  value  (note  that  7  =  —  <-»  /x  =  1).  Here  F  is  the  Gamma  function. 

3 

That  special  attention  is  needed  for  values  of  7  near  four-thirds  is  now  apparent  and, 

as  mentioned  earlier,  it  is  exactly  this  range  of  7  that  we  are  interested  in,  physically. 

Considering  c  = - a  to  be  small,  we  see  that 

4 

P(V  ~  2DoAl[(y^^)  -h  1  -f  0(e)]  (5.7) 

as  y  — +  0.  We  note  that  Pq^  (and,  in  fact,  Pq)  is  defined,  and  moreover,  continuous  as 
e  —  0.  It  is  convenient  to  recast  the  last  result  in  terms  of  z  to  ease  matching  to  the  next 
asymptotic  region  to  be  discussed.  The  result  is  that 

Po.  ~  2DoAlk^\k^{^^^)  +  -h  1]  (45) 

i-  4e  4f  J 


as  2  -+  0. 

The  solution  for  P  =  Po(y)  H - ,  will  continue  to  be  the  leading  order  behaviour  of  P 

until  contributions  from  the  ‘middle’-term  (proportional  to  P, )  become  leading  order.  The 
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location  of  this  next  region  of  interest  can  be  easily  located,  by  analogy  with  the  respective 
analysis  of  Section  3,  as  we  expect  c  >>  1  as  fc  — »  0.  Here  ^  ~  using  the  small-2 

asymptote  of  G,  and  we  expand 


P  =  2Do[1  +  k^6l^^P{g)  -h  . . .] ,  c  =  c62  +  •  •  • 


(5.9a, 6) 


where  62,  =  <^2(^1  f),  is  a  small  parameter.  The  P  form  is  implied  from  the  Pq  solution  for 
small  y.  The  new,  scaled  normal  variable,  y  is  taken  to  be  0(1)  and  is  defined  by 


2  =  <^2  (<<  1). 


(5.10) 


Substituting  into  the  pressure-amplitude  equation,  (4.37),  yields  the  equation  satisfied  by 


Pyy  -I- 


8go  p  ^  £0 
y{90  +  cy^) 


(5.11) 


where  yo  arises  from  the  small-2  asymptote  of  G.  This  can  easily  be  integrated  to  give 


Py{g)  =  {c  + 


fAf 


(C|/^  +  90)^ 


+  Dl. 


(5.12) 


with  Di  an  arbitrary  constant  to  be  determined  by  matching  with  the  previous  region.  As 
5  — »  00  we  see  that  Py  —>  Dic^,  thus  P^  D\^62^  and  matching  with  (5.8)  gives 


fcl  -  2o  _  1 


at  leading  order.  We  choose  D\  =  —Aqc~^  and  this  determines  the  size  of  c  for  small  k 


62  -  - 


(;S;T^)-4.  _  I 


>■ 

—k  >  +  *• 


(5.13) 


Taking  the  limiting  value  as  e  — ♦  0,  we  have  that  when  7  =  - 

— 1  >>  1 


(5.13) 


recalling  that  In  i  — »  —00  as  x  — »  0;  the  choice  of  I?i  <  0  ensuring  that  62  is  real  and  we 
allow  c  to  be  complex.  When  e  is  small  we  see  that  62  is  large  and  has  the  form  suggested 
at  the  start  of  this  section.  Thus  62,  as  given  by  (5.13),  leads  us  to  conclude  that,  for 
the  physically  realistic  values  of  7  around  four-thirds,  the  scaled,  relative  wavespeed  c 


increases  as  the  scaled  wavenumber  k  decreases,  this  growth  having  logarithmic  magnitude 
in  contrast  to  the  much  faster  algebraic  growth  found  in  the  no-shock  analysis  presented 
in  Section  3.  We  note  that  the  corresponding  frequency,  proportional  to  cfc,  tends  to  zero 
as  k  decreases.  The  k—  dependence  of  c,  for  fc  — >  0,  can  be  deduced  in  a  similar  fashion 
for  other  values  of  7,  and  we  note  that,  mathematically,  the  no-shock  result  of  Section  3 
is  the  limiting  case  7  — ►  00. 

All  that  remains  now  is  to  determine  c  and  this  is  easily  deduced  by  introducing  a 
further  region  where  z  ~  <<  1.  The  details  are  identical  to  those  of  Region  IV  in 

Section  3;  again  we  find  we  need  — »  0  as  y  — >  0,  resulting  in  an  expression  for  c 

r  sidy  .  Aj  „ 

K  +  flo)^ 


which  can  be  solved  for  complex  c;  after  an  elementary  contour  integration  and  a  little 
manipulation  we  obtain 


(5.14) 


which  corresponds  to  to  an  unstable  mode  of  (4.37). 

Thus  we  have  obtained  the  structure  of  the  vorticity  mode  in  the  strong  interaction 
region  at  small  values  of  the  wavenumber.  We  postpone  any  further  discussion  of  our 
results  until  the  following  section. 
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6  Discussion 

We  have  investigated  the  instability  of  flat  plate  hypersonic  boundary  layers  to  the 
vorticity  mode  of  instability.  This  inviscid  mode  is  associated  with  the  generalized  inflec¬ 
tion  point  of  the  basic  flow  and  is  thought  to  be  the  most  dangerous  mode  of  instability 
of  a  high  Mach  number  flow.  When  the  mode  is  neutral  the  wave  propagates  downstream 
with  the  speed  of  the  fluid  at  the  generalized  inflection  point.  At  wavenumbers  smaller 
than  the  neutral  value  the  mode  is  unstable  and  the  growth  rate  attains  it’s  maximum 
value  at  a  finite  value  of  the  wavenumber.  In  the  small  wavenumber  limit  the  growth  rate 
approaches  zero  and  for  the  non-interactive  boundary  layer  at  sufficiently  small  wavenum¬ 
ber  the  vorticity  mode  spreads  out  towards  the  lower  boundary  and  reduces  to  an  acoustic 
mode  at  a  countable  infinite  set  of  wavenumbers.  We  believe  that  a  similar  process  hap¬ 
pens  in  the  strong  interaction  case  since  there  the  acoustic  mode  is  correctly  described  by 
a  quasi-parallel  theory  there.  We  did  not  pursue  that  calculation  here  because  it  would  be 
essentially  unchanged  from  that  of  Section  3  except  that  it  would  be  made  somewhat  more 
complicated  by  the  necessity  of  treating  the  case  7  =  4/3  as  a  special  case  in  the  strong 
interaction  zone. 

We  believe  that  the  results  we  have  presented  in  Section  4  are  the  first  which  show 
the  effect  of  a  leading  edge  shock  on  any  form  of  hydrodynamic  instability.  Interestingly 
enough  the  shock  does  not  have  a  direct  influence  on  the  vorticity  mode;  thus  the  main 
effect  of  the  shock  is  to  restructure  the  boundary  layer  in  the  leading  edge  region  and 
thereby  influence  the  susceptibility  of  the  flow  to  inviscid  disturbances. 

The  vorticity  mode  eigenvalue  problem  was  formulated  in  the  interactive  region  along 
the  plate  at  0(1)  values  of  x.  However  to  the  authors’  knowledge  the  basic  flow  in  this 
regime  has  not  yet  been  calculated;  the  numerical  problem  was  set  up  by  Bush  (1966)  but  is 
sufficiently  difficult  to  have  remained  unsolved.  Thus  we  were  unable  to  solve  the  eigenvalue 
problem  in  this  regime  and  therefore  choose  to  consider  the  strong  interaction  regime  where 
a  similarity  solution  for  the  basic  state  is  available.  An  alternative  to  that  limit  would  have 
been  to  consider  the  weak-interaction  problem,  Bush  and  Cross  (1967),  where  a  different 
similarity  structure  holds.  We  choose  to  concentrate  on  the  strong  interaction  limit  because 
the  growth  rates  there  are  bigger  and  if  the  flow  is  indeed  unstable  there  the  stability  of 
the  flow  further  downstream  is  possibly  of  less  relevance. 

Unfortunately  we  are  unaware  of  any  experimental  observations  or  other  theoretical 
work  which  we  could  compare  with  our  results  for  the  strong  interaction  regime. 
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In  Section  3  we  showed  how  the  acoustic  inviscid  mode  emerges  from  the  small 
wavenumber  description  of  the  vorticity  mode.  Again  it  is  not  possible  for  us  to  com¬ 
pare  our  work  with  that  of  previous  authors  since  it  appears  that  the  finite  Mach  number 
calculations  available,  mostly  due  to  Mack,  have  either  being  carried  out  using  a  Chap¬ 
man  viscosity  law  or  a  combination  of  Sutherland’s  and  Chapman’s  law.  In  fact  Mack’s 
calculations  were  carried  out  using  a  combination  of  the  different  laws  so  as  to  efficiently 
model  the  viscosity-temperature  structure  of  the  fluid.  The  fact  that  the  calculations  of 
Cowley  and  Hall(1988)  and  Smith  and  Brown  (1989)  agree  so  well  with  Mack’s  calculations 
suggests  that  over  the  part  of  the  flow  where  instability  took  place  Chapman’s  law  was 
being  used;  in  the  case  of  the  vorticity  mode  this  is  clearly  a  bad  approximation  because 
the  mode  locates  itself  in  the  layer  where  the  basic  temperature  field  varies  rapidly. 
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solution  of  the  transition  layer  equation  (2.10a)  for  the 
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Figure  (3.6a)  The  real  part  of  the  wavespeed  Cq  as  a  function  of  L  for 
100  <  L  <  500  and  the  case  of  an  insulated  wall. 
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Figure  (3.6b)  The  imaginary  part  of  the  wavespeed  cq  as  a  function  of  L 
for  100  <  L  <  500  and  the  case  of  an  insulated  wall 
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Final  Report 

16.  Alien ect  The  loscabllity  of  hypersonic  boundary-layer  flows  over  flat  plates  is  considered.  The  viscos¬ 

ity  of  Che  fluid  is  taken  Co  be  governed  by  Sutherland’s  law,  which  gives  a  eucli  eore  accurate  repre¬ 
sentation  of  the  teaperature  dependence  of  fluid  viscosity  at  hypersonic  speeds  than  Cliapean's  ap- 
proxlute  linear  law;  although  at  lower  speeds  the  teaperature  variation  of  the  Bean  state  is  less 
pronounced  so  that  Che  Chapaan  law  can  be  used  with  soae  confidence.  Attention  is  focussed  on  Che 
so-called  "vortlclcy"  node  of  instability  of  the  viscous  hypersonic  boundary  layer.  This  is  thought 
to  be  Che  fastest  growing  inviscid  disturbance  at  hypersonic  speeds;  it  is  also  believed  to  have  an 
asyaptotlcally  larger  growth  race  Chao  any  viscous  or  centrifugal  instability.  As  a  starting  point 
we  investigate  Che  Instability  of  the  hypersonic  boundary  layer  which  exists  far  downscreaa  froa  the 
leading  edge  of  the  plate.  In  this  regiae  the  shock  chat  is  attached  to  the  leading  edge  of  the 
plate  plays  no  role,  so  that  the  basic  boundary  layer  is  non-lnceraccive.  It  is  shown  that  the  vor- 
ticicy  node  of  instability  of  this  flow  operates  on  a  significantly  different  lengchsoale  chan  that 
obtained  if  a  Chapaan  viscosity  law  is  assuaed  (see  Saith  and  Brown,  19B9).  In  particular,  we  find 
that  Che  growth  rate  predicted  by  a  linear  viscosity  law  overestlaates  the  size  of  the  growth  race  by 
O(M^).  Next,  the  developaent  of  the  vorcicity  aode  as  the  wavenuaber  decreases  is  described,  and 
it  is  shown  tliac  acoustic  aodes  eaerge  when  the  wavenuaber  has  decreased  froa  it's  0(1)  initial  value 
to  0(N-<).  Finally,  the  Inviscid  Instability  of  the  boundary  layer  near  Che  leading  edge  in  Che  in¬ 
teraction  zone  is  discussed  and  particular  attention  is  focussed  on  Che  strong  interaction  region 
which  occurs  sufficiently  close  to  the  leading  edge.  We  find  that  the  vorticlty  aode  in  this  regiae 
is  again  unstable,  and  that  it  is  concentrated  in  Che  transition  layer  at  the  edge  of  Che  boundary 
layer  where  the  teaperature  adjusts  froa  its  large,  0(m2},  value  in  the  viscous  boundary  layer,  to 
its  0(1)  free  streaa  value.  The  existence  of  Che  shock  indirectly,  but  significantly.  Influences 
the  instability  problea  by  apdlfylng  Che  basic  flow  structure  in  this  layer. 
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